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We consider long-range correlated disorder and mutual interacting particles according to a dipole- 
dipole coupling as modifications to the one-dimensional Anderson model. Technically we rely on 
the (numerical) exact diagonalization of the system's Hamilitonian. From the perspective of different 
localization measures we confirm and extend the picture of the emergence of delocalized states with 
increasing correlations. Besides these studies a definition for multi-particle localization is proposed. 
In the case of two interacting bosons we observe a sensitivity of localization with respect to the range 
of the particle-particle interaction and insensitivity to the coupling's sign, which should stimulate 
new theoretical approaches and experimental investigations with e.g. dipolar cold quantum gases. 
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I. INTRODUCTION & MOTIVATION 

When P. W. Anderson introduced a simple quantum 
model to represent a disordered lattice^ it turned out that 
the contained physics is surprisingly complex, i.e. there 
exists the phenomenon of Anderson localization which 
is related to an exponential decay of the quantum me- 
chanical probability distribution. A specific subclass of 
Hamiltonians Anderson studied reads 

L L-l 

Hs = ^ <^icjci + c+Ci+i + h.c. , (1) 

i=l 4=1 

which represents a single particle in a one-dimensional 
chain of L sites with random (onsite) potential e; and ki- 
netic energy J {hopping energy). We assume the formalism 
of second quantization where the operator c|^^ annihi- 
lates (creates) a boson at site i and hence we refer to Eq. 
(1) as the disordered Bose-Hubbard model (without mu- 
tual interaction). One Can either derive Hs from the tight 
binding approximation of a continuous modeP of non- 
interacting particles in an external potential or one ab 
initio takes it as a discrete model. 

Based on the renormalization group flow"'^ idea one 
can argue that the conductance of a disordered solid may 
vanish for sufficient large systems^ . In particular the one- 
dimensional disordered system, Eq. (1), becomes an in- 
sulator in the limit i — >■ cxd. 

An explicit argument for localization of all states \E) 
satisfying Hs\E) = E\E) can be established by ex- 
ploiting the transfer matrix method^ and a theorem due 
to Fiirstenberg'' . On the other hand there is Bloch's 
theorem*^ which induces periodic (delocalized) states for 
a periodic lattice potential e^. Hence the random nature 
of the potential must be obviously the key feature that 
leads to localization: A criterion based on the differen- 
tiability of the disorder potential was recently studied 
to understand the degree of randomness necessary for 



delocalization^. 

In order to investigate the impact of correlation on 
localization we introduce a specific disorder model that 
extrapolates from a pure random sequence ei, e2, . . . , cl 
to a periodic, and thus correlated, structure in Sect. II. 
It was first used by Moura & Lyra® and similar investiga- 
tions followed^ In a first step we will review the model 
on the basis of three different localization measures. Fur- 
thermore, we utilize one of those quantities to establish 
the phase diagram — namely, the dependence of local- 
ization with respect to disorder strength and the amount 
of correlation among the e^. While a previous study by 
Shima et al.^^ focused on the properties of states in the 
band center, our measure in use accounts for the global 
aspect, i.e. it incorporates properties of the entire spec- 
trum. 

The potential interest in the sensitivity of Anderson 
localization on correlated disorder arose due to recent ex- 
periments with Bose-Einstein condensates (BECs) where 
the direct observation of the atomic density distribution 
provides access to the quantum mechanical probability 
distribution^^'^"^. Besides correlation the impact of inter- 
action is an issue one is naturally faced with when study- 
ing localization in BECs. Anderson already mentioned 
the importance of particle interaction^"' on localization 
and worked out theoretical investigations in collabora- 
tion with L. Fleishman '^'''. Over the years various aspects 
and features of interacting particles in a random potential 
have been figured out, but the problem remains a chal- 
lenging topic for present research since results from dif- 
ferent approaches do not always coincide. Just recently, 
the phase diagram of the 3-dimensional disordered Bose- 
Hubbard model was established^^. 

In our discussion on interacting particles in the pres- 
ence of a disordered onsite potential (Sect. Ill) we first 
want to focus on a suitable definition of multi-particle 
localization and then turn to the extension of Eq. (1), 
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namely 



s 



(2) 



with a (two-body) interaction potential f/y- . We model in- 
teraction according to a magnetic dipole-dipole coupling 
obtained in BEC experiments with dipolar gases^''^^ 
which is beyond the standard treatment of the onsite 
interaction term Uq ~ 1) (with fii = cfci) present 

in the disordered Bose-Hubbard model. By explicitely 
diagonalizing Hmp for two interacting bosons we explore 
the relevance of interaction (especially its range) for local- 
ization and discuss an interesting symmetry involving the 
sign of the interaction potential Uij. 



II. CORRELATED DISORDER 

To model correlated disorder we use the following pre- 
scription for the onsite potential values in Eq. (1): 
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27r , , \ ( a>0 
iV^^ + '^V 1 N~»L ■ 

(3) 

Here, a denotes the correlation parameter, (pk G [0, 27r) 
are ^ uniformly distributed random phases and N is 
a natural number that should be much larger than the 
number of lattice sites'**^ L. Roughly speaking, the 
(i = 1 . . .L) represent the discrete Fourier transform of 
k~°', i.e. an algebraic decaying power spectrum^^ . One is 
therefore used to refer to the ei as long-range correlated 
and indeed if we turn back to the continuum limit, we can 
argue that the are correlated according to an algebraic 
decay for a G (0,1). The case a — Q corresponds to 
almost uncorrelated disorder which is close to the perfect 
disorder Anderson assumed in his model. More details 
on properties of the from Eq. (3) are provided in App. 
Al. 



A. Localization Measures 

In order to detect the localization-delocalization transi- 
tion we introduce three different measures: 

a) the Normalized Standard Deviation^'^ (NSD) 

b) the Inverse Participation Ratio^^ (IPR) and 

c) the Nearest Neighbor Distribution^^ (NND). 

While the NSD and the IPR are derived from the spa- 
tial/site probability distribution 



'il?E{i) = \{i\E)\^wit\\Hs\E)^E\E) and 



|0) 



(4) 

where Ci |0) =0, the NND depends on the spectral values 
E only. Given a chain of L sites we calculate 



where (.) = {E\.\E) denotes the expectation value and 
hence the NSD derives the (spatial) variance of the lattice 
site index i for some \E) with respect to the state 
whose probability distribution is uniform on the lattice. 

In the case of the IPR one quantifies the (inverse) num- 
ber of sites where significantly differs from zero. 
Since X^i^i i^si'^) = 1 we state that 



IPRiEj^Y.i^'M'' 



(6) 



with L defined as the number of sites that are occupied 
by \E). Indeed, we can convince ourselves that IPR ^ 1 
and IPR ~ l~+oo^ ^ localized and delocalized 

states, respectively. 

In contrast, the NND considers the spectral properties 
of Hs- More precisely, we evaluate fluctuations of level 
spacings around a local mean s„ by computing 



{En+1-En)/Sn with S„ 



E^, 



En- 1 



TO + 1 



(7) 



and deriving the distribution P{s) of the nearest neigh- 
bor spacings Sn- Here, we labeled the spectral values 
according to Ei < E2 < ■ ■ ■ < El and the division by 
Sn unfolds'^ the level spacings to relate different parts of 
the whole spectrum to each other. The procedure of un- 
folding is not unique^'' and to is left as a free parameter 
that defines the notion of local. We choose it such that 
TO ^ L on the one hand and to include enough energy 
values for reasonable statistics on the other hand; in fact, 
we used to ^ 20 for L ^ 10"^ . 

The NND is not as obvious as the former measures^'^'^'"'. 
Intuitively, the argument works as follows: Taking two 
different energy eigenstates \E) , \E') that do not signif- 
icantly overlap (in site space) we assume them to be al- 
most orthogonal, i.e. they are in some sense independent 
from each other (localized) and nothing prevents the corre- 
sponding energy values to be arbitrary close (s -s- 0). But 
if the overlap increases (extended states) the levels start to 
repell"^^: E ^ E' . Quantities that are statistically in- 
dependent exhibit a Poissonian distribution and thus the 
corresponding NND should be'^'' P{s) ^ e~*. In contrast, 
P(0) = is reasonable to expect for delocalized states. 



B. Numerical Results 

Turning back to Eq. (1) we rescale Hg by J, i.e. Hg — > 
J~^Hs, which does not alter the Hamiltonian's eigen- 
states but multiplies the spectral values E by a. factor 
of J~^. Since we restrict the Ci to the (finite) interval 
[—y, y] we introduce the parameter 



A/J 



(8) 



(5) 



that indicates the strength of disorder. Thus we are faced 
with a two-dimensional set of system parameters (a, k) 
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LL at K = 1 
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FIG. 1: (color online) Energy-resolved plot of the localiza- 
tion length (LL) for different correlation parameter values a 
to compare our results from exact (numerical) diagonalization 
to investigations that use the transfer matrix method to ob- 
tain the Lyapunov exponent 7 (cf. Ref. [8], fig. 5). We used a 
system of L = 3500 sites and averaged over 50 disorder real- 
izations. The inset (transparent to the main panel in the back- 
ground) presents numerical data of LL(Q) for different system 
sizes which underpins the reported delocalization transition 
around a — 2. 
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FIG. 2: (color online) Plot of the distribution P(NSD) of 
the (spatially) normalized standard deviation Eq. (5) at 
intermediate disorder k — 1 (J = A) to support that 
localization— ^-delocalization takes place rather smoothly than 
according to a sharp transition. As for Fig. 1 we used 
L — 3500 and the averaging was taken over 50 disorder real- 
izations. The inset (transparent to its background) shows results 
for different system sizes and supports our crossover picture 
as well (details in the main text). 



and, from the following as well as App. Al, it becomes 
clear that increasing a corresponds to increasing corre- 
lations up to long-range correlated disorder. Moreover, 
since the NSD and the IPR only depend on \E), we 
rescale and shift the spectrum such that the rescaled val- 
ues E ^ e obey —0.5 < e < 0.5 when plotting these 
quantities resolved in energy. 

To set the stage, we want to relate our numerical re- 
sults to previously published ones. Since it is common 
practise to attempt localization by computing the Lya- 
punov exponent^^ we establish an exponential fit to 
\E). More precisely, our numerics picks out the maxi- 
mum -01; (iq) and performs an exponential fit into the di- 
rection singled out by Taax.{if),L — io). From ip'^ /it(*) ^ 
exp (—7 |i — iqI) we extract 7(i?) and define the localiza- 
tion length 



LUE)^^-\E)/L 



(9) 



The corresponding result for intermediate disorder (k = 1) 
is shown in Fig. 1 and it exhibits reasonable qualitative 
agreement with Ref. [8], fig. 5; namely, delocalized states 
(7~i~L->-LL~l) arise in a finite range around the band 
center e = when the correlation within the disorder po- 
tential Ei is increased. This interpretation also coincides 
with the plot presented in Ref. [9], fig. 3, where localiza- 
tion was quantified by means of a measure based on the 
density of states 



DOS(£;) 



dn 
dE 



(10) 



which defines the number dn of energy eigenstates \E) in 
a given interval [E,E -\- dEf^^ . 

Moreover, the work of Moura & Lyra mentioned above 
supports delocalization at a = 2 for states at the band 
center. Therefore we depicted a cut of Fig. 1 at e = 
in the relevant correlation parameter range a G [1,3] 
and checked LL(0) according to a finite size analysis, in- 
set of Fig. 1. Indeed, up to a « 2 the LL(0) decreases 
with increasing system size L while above this correlation 
parameter value it remains relatively constant. Assum- 
ing that this trend persits for even larger system sizes, 
the data support that a = 2 marks a qualitative differ- 
ence between systems with smaller and larger correla- 
tion, respectively: While the relative extend of the states 
with respect to the system size falls off (localization) for 
< ^ 2, it remains constant (extended states) for a > 2 
in the thermodynamic limit L — ^ 00. 

However, we would like to address the question where 
(in terms of a) does the localization-delocalization tran- 
sition takes place from the perspective of the measures 
introduced in Sect. II A: Instead of a sharp transition 
at a = 2, which is supported by Fig. 1 and publica- 
tions mentioned above, we suggest a smooth crossover 
in 1 < a < 2. Investigating the distribution P(NSD) 
most obviously illustrates this statement and we present 
it in Fig. 2. Although the major fraction of localized 
states (NSD < 1) becomes delocalized at a « 2 there 
is a finite fraction that splits apart the localized region 
(0 < NSD < 0.2) around a = 0.5 and drifts towards the 
delocalized regime (NSD as 1) up to a w 2. This observa- 
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FIG. 3: (color online) Smooth localization-delocalization 
crossover from the perspective of the simultaneous analysis 
of the IPR and the DOS. Increasing correlation attracts the 
system's states from the band edges (|e| = 0.5) to the band 
center (e = 0) where delocalized states are present and thus 
the system becomes more and more delocalized. Finally, lo- 
calized states almost completely vanish above a = 2. The 
inset of the upper panel (transparent to the background again) 
establishes the finite size analysis similar to the inset of Fig. 
2. 

tion remains stable in the numerically studied range of 
system sizes (L = 1000 . . . 3500). 

An explicit plot to the finite size analysis performed 
with respect to our crossover picture is given by the inset 
of Fig. 2. Convoluting the NSD with the DOS, i.e. 



NSD = J dE NSD(i;)DOS(£;) (11) 

yields a single quantity to characterize the global local- 
ization property of the disorderd system at fixed correla- 
tion parameter a. In the spirit of the previously applied 
LL-dependence on L, the inset of Fig. 2 provides the sys- 
tem size dependence of the averaged normalized standard 
deviation NSD: Up to a w 1, NSD decreases with in- 
creasing system size L indicating that the Hamiltonian's 
eigenstates stay localized in the thermodynamic limit. 
Above a ss 1 the averaged standard deviation that is 
normalized to L remains unaltered for the accessed num- 
ber of sites — a signature of delocalization. Remarkably, 
the NSD forms a plateau in 1 < a < 2, before growing up 
in absolute value with increasing correlation. One more 
evidence to the crossover picture we are in favor of. 
Similarly, we investigate the IPR taking into account 



the DOS from Eq. (10). Since we are interested in the 
localization property of the whole system it is necessary 
to account for the number of states with a certain local- 
ization measure value. Hence the IPR and the DOS are 
plotted simultaneously in Fig. 3 and the inset of the up- 
per panel presents an averaged version of the inverse par- 
ticipation ratio in total analogy to Ec^. (11) with NSD(iJ) 
replaced by IPR(£'). 

From a naive point of view we could conclude that 
even for weak correlation {a < 0.5) the system is delo- 
calized since there are only a few localized states near 
the band edges (|e| = 0.5) while the rest of the spectrum 
exhibits an IPR value corresponding to delocalization. 
Moreover, increasing correlation enhances localization up 
to a « 1 which seems unexpected. But if we include 
the DOS (lower panel of Fig. 3) the picture drawn before 
reveals: We realize that almost all states have energies 
near to the band edges for a < 0.5 and thus the system 
is localized. The crucial fact leading to the localization- 
delocalization crossover is that for increasing correlation 
the band center attracts the states from the band edges 
to the regime where delocalized states are permanently 
present (|e| < 0.3). Localized states at the band edges 
follow that trend only up to a « 1 and for stronger cor- 
relation they start to disappear until they are almost van- 
ished at aw 2. According to an argument presented in 
App. A 1 it is reasonable to refer to disorder with a > 1 
as highly correlated and therefore it provides a hint to 
the suggestive importance of a = 1 in our data. 

Again, this observation is confirmed by a finite size 
analysis as shown in the inset of the upper panel of Fig. 
3. We note that the IPR-values stay quite unaltered for 
a < 1 when L is increased: Since IPR ~ ~ 
(cf. Eq. 6) corresponds to extended states which tends to 
zero in the limit L — >■ c», the decreasing behavior of 
IPR(i) for a > 1 suggests the existence of delocalized 
states. The bump in the crossover interval a £ [1,2] 
may represent a finite size effect that eventually vanishes 
in the thermodynamic limit. 

Finally, we would like to take a closer look at the 
smooth crossover from the perspective of the energy spec- 
trum. Therefore we investigate the NND from Sect. II A 
which also prepares the relevant measure we will use to 
quantify the {k, a)-dependence of localization. Plots of 
the NND for different correlation parameters a are pre- 
sented in Fig. 4. For small a the result is near to the 
Poissonian distribution e"**. Deviations from the theo- 
retical prediction are due to the finite system size L and 
the fact that a = does not exactly correspond to uncor- 
related but weakly correlated disorder (see App. Al). We 
further note that the data points for a — and a — 0.5 
almost coincide, but when increasing a above 0.5 the dis- 
tribution starts to deviate more significantly from e"'' to 
become peaked around 5 = 1. Returning to Eq. (3) we 
note that for a 3> 1 the summation over the (discrete) 
momenta k picks up less cosine modes which leads to a 
periodic onsite potential in the limit a — oo. Hence the 
spectrum of becomes more regular for increasing a 
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FIG. 4: (color online) Semi-logarithmic plot of the nearest 
neighbor spacing distributions for several correlation parame- 
ter values Q (L = 3 ■ 10*, k = 1). The broken line is to compare 
them to the Poissonian distribution Pp which we expect for 
localized states. As apparent, the distribution does not con- 
verge to PwD of the Gaussian orthogonal ensemble'^ since for 
a ^ 1 the onsite potential becomes periodic, local fluctu- 
ations among the Ei vanish and thus P{s) must be peaked 
around s = 1 (s„ = (i?„+i — En)/sn 1, a > 1). A more 
detailed discussion is given in the main text. 



vanishes except for its diagonal elements Hs^a = ti and 
the first off-diagonal entries Hs^u+i = Hg^i^u = J when 
open boundary conditions (cf. App. A 2) are imposed. In 
the limit where the (uncorrelated) disorder strongly dom- 
inates the kinetic energy (A > J), the single particle 
Hamiltonian becomes approximately diagonal to zeroth 
order in ^ and the solution to the Schrodinger equation 
reads \Ei) sa \i) with Ei w e^. Hence the almost purely 
random/uncorrelated sequence of eigenvalues Ei should 
obey a Poissonian distribution in the NND. 

However, there is a crucial fact we would like to keep in 
mind: Although the Hamilton matrix of is real valued 
and symmetric (in site space) — since we deal with static 
disorder and real valued constant hopping energy J — 
the elements Hg ij do not fulfill one of the conditions'"^ 
from which PwD is derived as the corresponding NND of 
quantum chaotic systems. Therefore we will not natu- 
rally expect the Wigner surmise to arise in the presence 
of delocalized states. Indeed, Fig. 4 numerically verifies 
this conjecture. 

An issue closely related to the analysis of the NND is 
the question on a critical distribution Pc that marks the 
transition from localization to delocalization. One way 
to detect this transition by means of the NND is to argue 
for a suitable intermediate distribution between Pp and 
PwD as done in e.g. Ref. [27]. Usually one investigates 
the quantity 



and (local) fluctutations eventually vanish for sufficiently 
strong correlation due to Bloch's theorem. From the def- 
inition of the unfolded level spacings s„, Eq. (7), we con- 
clude that all s values should be centered around 1 due 
to the normalization by the local mean s„. 

At this stage it is maybe advisable to discuss the 
Wigner surmise PwD — a specific NND connected to 
random matrix theory'^ — in the context of localization, 
especially in our specific model system under considera- 
tion. In the literature^^"^^ it seems to be common prac- 
tise to associate 

Pp(s) EEexp(-s) (12) 

with localized states and signatures of level repulsion 
{P{s < 1) < 1) in the NND with extended ones after 
having unfolded the spectrum of the corresponding dis- 
ordered system. In the case of a time reversal symmetric 
quantum chaotic system — represented by real, symmet- 
ric random matrices — the NND adopts the Wigner sur- 
mise 

Pwd(s) = |sexp(-^s2^ (13) 

in the thermodynamic limit (cf. Ref. 4 as well as references 
therein) . 

As we scetched it in Sect. II A it is quite plausible to 
adopt the first correspondence. Indeed, the matrix repre- 
sentation oi Hs in the site basis {]«)} becomes tridiagonal, 
i.e. 

(14) 



_ Jds [P(,s)-Pwd(s)] , . 

/ds [Pp(.)-Pwd(s)] ^ ' 

based on the NND P(s) where the integrals may be cho- 
sen (including a proper argumentation) On a Suitable interval 

where one expects sensitivity of P(s) on localization. Ac- 
cording to our discussion from Sect. II A one may take 
s € [0, Smax] with Srnax Smaller than the smallest root 
of Pwd(s) — Pp(s). But since our one-dimensional dis- 
ordered system in use does not converge to the Wigner 
surmise it is rather vague to apply Eq. (15) in order to 
resolve the question on the impact of correlation on lo- 
calization. 

Nevertheless, the previous discussion invites us to in- 
troduce a similar measure for plotting the dependence 
of localization on the parameter space (a, k) which will 
be the last purpose of this section on correlated disor- 
der. We determine the deviation of P{s) from Pp as a 
quantity for the degree of delocalization. Thus we define 

^nndI'*:") = SmL / ds[P(s,K,a) - Pp(s)]2 (16) 
Jo 

where Smax introduces the numerically necessary restric- 
ti on to a finite range of s-values and we use ^nnd = 
\/ -^NND ^ ^ localization measure. Numerically, Xnnd 
features the striking advantage that it is sufficient to just 
compute the eigenvalues of the Hamilton matrix Hs^ij 
and therefore we are able to reach larger system sizes L 
with our code. 



6 




disorder strength k 



FIG. 5: (color online) (k, a)-dependence of localization in the 
one-dimensional Anderson model from the perspective of a 
measure based on the nearest neighbor spacing distribution 
(NND) of the Hamiltonian's spectrum with L = 3 • 10*. The 
finite size analysis with the measure A[Lnnd]^: defined in Eq. 
(17), is shown as an inset (white). Furthermore we included 
the region which corresponds to extended states according to 
the work of Shima et al.^^ who tackeled the phase diagram by 
means of states in the band center E — using the transfer 
matrix method. The (•) mark the minimal deviation of the 
NND from Psp (cf. Eq. (18)) for fixed disorder strength k and 
the (x) indicate the same minima with respect to Pwd . 

As we know from the previous considerations on the 
locahzation-delocahzation crossover for intermediate dis- 
order (k = 1, first part of Sect. II B), correlation with 1 < 
a < 2 seems to be important for delocalization; vary- 
ing a in [0, 2.5] will cover it. In the limit k — ^ 0, i.e. 
vanishing disorder, it is obviously hardly possible to re- 
cover Anderson's result for infinitely extended systems 
and arbitrary small disorder, since we have to face the 
fact to be generally restricted to finite system sizes L 
and finite numerical precision. Therefore we investigate 
the domain (a, k) G [0, 2.5] x [0.5, 10] which spans a wide 
range of disorder strengths k for the interesting amount 
of correlation among the onsite disorder potential val- 
ues ei. The corresponding result is shown in Fig. 5 with 
L = 3 ■ 10* and the inset presents the finite size analysis 
we performed: Let AL denote the (fixed) difference be- 
tween two system sizes L and L' such that L = L' + AL 
and we define 

^{Lnnd?{L,AL) = \Ar^ [ dnda x 

J A 

[L'^j^D{K,a,L) - L'^j^^{K,a,L')]'^ (17) 



with 1^1 = I [0,2.5] X [0.5, 10] I = AkAo the area of in- 
tegration. The white inset plots that quantity versus L 
having fixed AL = 1000; it becomes clear that Fig. 5 
exhibits a convergent trend for increasing L. 

Before turning to a more detailed analysis of 
iNND('*,ct) let us come back to our comments on a crit- 
ical distribution P^. As stated before we do not expect 
the NND to converge to Pwd in the limit a — cx) (delocal- 
ized system). In our case it is not well justified to assume, 
e.g., the Semi-Poisson law"^^ 

Psp(s)EE4sexp(-^2s) (18) 

that resembles the linear increase of Pwd(s) for s ^ 1 
and an exponential decay for s ^ 1 to be critical. Nu- 
merical evidence on our claim provides Fig. 5 where we 
plot both the minimal deviation of the NND from Pgp 
(•) and Pwd (x) for fixed disorder strength k. The notion 
deviation is defined in total analogy to Eq. (16). In the 
case of the convergence of the NND to Pwd we would 
expect the minima (x) to appear at the maximal a-value 
simulated — independent from the disorder strength k. 
But, as visible in Fig. 5, the NND passes by Pwd below 
a ~ 2.2. But surprisingly, the deviation minima of the 
NND from Psp (•) appear to lie in the crossover region 
1 < a < 2 we emphasized during our previous numerical 
studies. Moreover, they follow the qualitative trend of 
our localization measure Lnnd which one may naively 
exploit to argue in favor of Psp as a suitable critical dis- 
tribution. But a closer look on the precise shape of the 
NND with minimal deviation from Pgp reveals signifi- 
cant differences from the Semi-Poisson law that become 
even more obvious with increasing disorder strength k. 
A corresponding plot is shown in Fig. 6. Therefore, we 
do not support the Semi-Poisson law Pgp as a suitable 
critical distribution. In fact it is designed as a hybrid 
between the two limiting cases of a Poisson distribution 
Pp and the Wigner surmise Pwd which does not apply 
to extended states in our specific model. 

However, there has been a similar investigation to 
our phase diagram in Ref. [11] and we would like to 
discuss some conclusions one may draw when compar- 
ing the two results now. Shima et al. report numerical 
studies that support the existence of extended states for 
(k, a) G [0, 4] X [2, 5]^'^ on basis of a quantity A similar to 
our LL evaluated at the band center e = (cf. their Fig. 
6). We separated the corresponding region of extended 
states in Fig. 5 by a dashed line. There are two main dif- 
ferences we would like to point out with respect to our in- 
vestigations of Lnnd('*j et)' First of all, the localization- 
delocalization crossover depends on the disorder strength 
and second Lnnd varies over approximately two orders 
of magnitude in 1 < a < 2 quite independently from the 
value K. Therefore the evaluation of the NND underlines 
the picture of a smooth crossover between localization 
and delocalization drawn before. The qualitative result 
of delocalization for all /c-values under consideration is 
in accordance with our line of reasoning that the disor- 
dered system has to eventually delocalize for a — > oo due 
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FIG. 6: (color online) NNDs for different disorder strengths K (L = 3 ■ 10^) to illustrate that the Semi-Poisson law Pgp as a 
hybrid between the Poisson distribution Pp and the Wigner surmise Pwd is not an appropriate critical distribution to indicate 
delocalization for our model under consideration. The NND P{s) with Qmin (0) is singled out by its minimal deviation from 
Psp for fixed k as described in the main text (cf. also (•) in Fig. 5). Plotting P{s) with Cfmin ± 0.1 underlines the two statements 
we emphasized: a) decreasing values of the NND for s — >■ (level repulsion) and b) development of a peak of P{s) around s = 1 
for increasing correlation parameter a due extended states (increase of P(l)). The inset in the right panel shows the dependence 
of the minimal deviation from Psp with respect to the disorder strength. 



to Bloch's theorem. Therefore it would be interesting to 
push the numerical analysis of Shima et al. beyond their 
maximal correlation a -H- p = 5. However, assuming that 
the deviation of the results of Shima et al. and our in- 
vestigations is not caused by some hidden technicality'''*, 
the discrepancy of our conclusions to Shima et al. may 
reveal two useful lessons: 

1. It reminds us of the difhculty of a proper interpre- 
tation of the NND with respect to localization. The 
quantities A and LL, respectively, are much closer 
to the original picture of localization drawn by P. 
W. Anderson and thus one may prefer it in cases 
where physical intuition is hard to gain. The NND 
always needs reasonable understanding of the sys- 
tem under examiniation. 

2. While Shima et al. focused on the states in the band 
center of the disordered one-dimensional system, 
the NND indirectly takes into account the (averaged) 
localization property of all states. Assuming that 
both quantities properly account for localization, it 
seems that, at least for (k, a) G [4,8] x [2,5], the 
properties of states with E = are not dominating 
enough to show up when all states are included. 

As a closing remark we would like to note that the pre- 
vious discussion points out that the application of the 
NND as a localization measure is not as straightforward 
as the other quantities in the particular system we are in- 
vestigating. Even though its appealing property of being 
basis independent and thus easily employable for study- 
ing interacting particles, its interpretation remains com- 
plicated. Eventually, this observation prevents us from 
utilizing it in order to detect localization of multiple par- 
ticles and we would rather like to argue for an intuitive 
observable in the section below. 



III. LOCALIZATION IN THE PRESENCE OF 
INTERACTION 



In this second part of our discussion on delocalization in 
the one-dimensional Anderson model we want to turn to 
the question how localization is effected by the presence 
of interaction. This problem basically goes beyond the 
scope of the quantum physics of Hs originally encoun- 
tered by Anderson, but — as mentioned in the begin- 
ning — the question of the impact of interaction already 
attracted his attention in the late 1970s. Nevertheless 
the problem has remained a challenging topic and exper- 
iments with BECs in optical lattices have re-raised the 
focus on it during the past few years. 

With Eq. (2) at hand, namely 



L L-l L 

2—1 i— 1 — l 

(19) 

we introduce an interaction term Uijcf CiCj where for 
Uij — UoSij one encounters the disordered Bose-Hubbard 
model assuming bosonic commutation relations for the 
lattice site annihilation (creation) operators c^^K Uq = 
const, describes contact/ onsite interaction^ i.e. particles 
interact only when occupying the same site. But as we 
will carry out in Sect. IIIB it is convenient to assume 
the more general case, Eq. (19), to describe experiments 
with either Rydberg gases'^*''^^ or (dipolar) BECs*^ where 
particles interact even if separated by a (large) number of 
lattice sites. In particular, we use an algebraic decaying 
interaction potential and therefore conceptually bridge 
from long-range correlation to long-range interaction. 
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A. Defining Localization 

Before proceeding we have to deal with the question on 
how to treat locahzation for multiple (interacting) parti- 
cles, since localization was originally defined for the sin- 
gle particle problem. The literature provides several ap- 
proaches for a suitable definition. E.g. one can use the 
Hausdorjf measure'''^ as a distance between two states 
ipo and ipt where the latter is the former one evolved in 
time by the Schrodinger equation ihdt = Hmp {ip}- In 
analogy to Anderson's initiating paper^ one can then ask, 
roughly speaking, for the absence of diffusion in terms of 
that distance. Another idea is to relate localization to 
macroscopic observables and associate the notion of lo- 
calization to vanishing electrical conductivity'^'*. Further- 
more one can directly study localization in Fock space as 
performed, e.g., in Ref. [35]. An ansatz independent of 
the specific form of the Hamiltonian — and thus easily 
applicable to the interacting problem — is the analysis 
of spectral statistics (cf. Eq. (7))^'' which is driven by the 
analogy between random matrices in quantum chaos and 
the random Hamilton matrices due to disorder. But as 
extensively discussed in Sect. II B its strength of being 
representation independent may come along with some 
difficulties due to a proper interpretation of the corre- 
sponding results. Furthermore, the application of the 
NND to systems of multi-particles shifts the argument of 
level repulsion given in Sect. II A to localization in Fock 
space. 

For our purpose we want to propose a rather intuitive 
solution inspired by a quantity which seems to be natu- 
ral under the scope of experiments with cold atoms and 
which is connected to the particle density. Therefore we 
try to motivate a projection of the many-particle state 
to the lattice sites i = 1 . . . L and apply a one-particle lo- 
calization measure afterwards. Of course, the procedure 
has to coincide with the notion of single particle localiza- 
tion if we return to a set of Fock states with one particle 
at a given site. 

Since one is able to directly image the density profile of 



Bose-Einstein condensates in experiment 



consider 



the question on the probability pi to find at least one par- 
ticle at site i. Using the language of second quantization 
we define the projection operator 



the site basis \i) by the probability distribution pi{E). 
For the single particle problem this projection is ex- 
actly given by the site basis representation of |£'), i.e. 
Pi{E) = |(z|£')|^ whose properties we studied for corre- 
lated disorder above. Hence the natural step to perform 
next is to apply the NSD or the IPR on pi {i = 1 . . .L). 
But there is an argument excluding the application of 
the NSD measure: Suppose, two particles are localized 
at certain sites ii and ?2: i-e. Pi is peaked for i £ {ii,i2}- 
If |ii — 12! ~ L the NSD will yield values of order 1 and 
we would refer to the corresponding \E) as delocalized. 
Thus we exclusively consider the IPR and our definition 
of localization turns, roughly speaking, to the associa- 
tion on how many sites are occupied by the interacting 
particles. 



B. Modelling Interaction 

As announced above we want to consider long-range cor- 
relation inspired by exeriments with dipolar gases. Since 
no distance Ax between neighboring sites enters the the- 
ory we can think of the limit L — > 00 twofold. It can ei- 
ther refer to approaching an infinitely large system with 
finite Aa; or the continuum limit {Ax — 0) of a finite total 
system length x — LAx (like in experiments). The following 
discussion prefers the second picture. 
To properly model Uij let us consider 



U^j = U{1) ^ i- 



for / 



(22) 



according to the dipole-dipole coupling. However, for I 
1 the interaction has to be renormalized to reach a finite 
value at Z = which corresponds to the (Bose-Hubbard) 
contact interaction strength Uq- To respect these limits 
we use the L-dependent'^'' interaction 



± 



I 



(23) 

Note that A/ > specifies the range of interaction and 
the sign ± determines repulsion and attraction, respec- 
tively. 



(20) 



where is the eigenvalue of the number operator hi = 
cf Ci counting the number of bosons at site i and the 
Kronecker delta Ski ensures that Vi is well defined for 
Ui = 0. Given an eigenstate \E) of H^^p we can answer 
our question from above by 



p,{E)^\{E\iV+r^)\E)\ 



(21) 



and the constant of proportionality is obtained by the 
normalization condition 1 = ^f=iPi- Therefore we end 
up with a projection of the multi-particle state \E) onto 



Two bosons with dipole-dipole coupling: 
Numerical results 



So far our considerations to define localization and the 
type of interaction were rather general on the total par- 
ticle number n = X^iLi "-i ^-nd the system size L, respec- 
tively. But for numerical simulations we need to set up 
a practically tractable situation: Since the dimension of 
the corresponding Fock space grows exponentially with 
the total number of particles when assuming constant fill- 
ing"^^, we investigate the specific case of constant particle 
number n. Therefore the complexity class of dim Hmp 
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shrinks to 0{L'^) and in the specific case of two particles 
we obtain 

dimff^p^ + =0{L^) , n = 2 . (24) 

Hence we are able to reach much larger system sizes L 
in contrast to the situation L ^ n and we are especially 
interested in the (well-established) toy model'^^"''^ of two in- 
teracting particles since D. L. Shepelyansky and Y. Imry 
provide arguments for the invariance of localization with 
respect to the sign of interaction which we can check nu- 
merically. Of course, we do not expect that studying 
just two interacting particles in a disordered potential 
will fully account for the effects of finite densities in cold 
atom experiments, but before turning to full complex- 
ity we may gain intuition on the problem by means of 
this academic example. As we will describe below there 
is some phenomenological reasoning that relates delocal- 
ization by correlation and interaction. 

As in Sect. II B we apply exact numerical diagonaliza- 
tion to the Hamiltonian Hmp and hence we have to spec- 
ify a suitable basis: We take {\ni . . .rii . . . n^)} and App. 
A 2 provides details on the explicit form of the Hamilton 
matrix for arbitrary size L{L -\- 1)/2. The inclusion of an 
interaction term Uij to the problem enlarges our space of 
system parameters from (k, a) to (k, a, ztuo. A/), where 

uo = |i7o|/J (25) 

accounts for the rescaling discussed in the beginning of 
Sect. II B. We should concentrate on the new degrees 
of freedom (±uo,A/) and therefore we fix k as well as 
consider uncorrelated disorder a — Q where Anderson 
localization of all states is a proven fact in one dimen- 
sion for vanishing interaction. To obtain a measure 
that characterizes the localization property of the whole 
system (parameters k,o fixed), i.e. for all \E) satisfying 
Hmp\E) = E\E) (cf. Fig. 3) we convolute IPR(i;) ex- 
tracted from Pi{E) with DOS(i?) as performed several 
times during our numerical studies in Sect. II B (cf. Eq. 
(11)): 

TPR{±uo,Xi) = J dE IPR{E)J:)0S{E) . (26) 

The corresponding result is shown in Fig. 7. As inde- 
pendently predicted by Shepelyansky and Imry local- 
ization seems to be invariant under the transformation 
Uij — ?> —Uij. Moreover, the range A/ determines whether 
the increasing onsite interaction strength uq (weakly) de- 
localizes or tends to localize the two bosons. This effect 
sets in when |J7o| ^ J, i.e. when the interaction energy 
becomes comparable to the kinetic contributions in Hmp- 
A plausible argument for the observation that our re- 
sults yield independence from the sign of the mutual 
particle interaction relies on a discrete (spatial) symme- 
try of the Hubbard Hamiltonian''^^'^"'^ which we apply to 
Eq. (19). Suppose we transform the spatial wave func- 
tion {j\E) such that all, say, odd site j contributions are 
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FIG. 7: (color online) DOS convoluted inverse participation ra- 
tio IPR on the basis of the projection, Eq. (21). The quantity 
was evaluated for several (rcscaled) onsite interaction strengths 
uq and interaction ranges A/ at fixed k = 5. It is remarkable 
that the result seems to be independent of whether the dipole- 
dipole coupling is repulsive (orange) or attractive (blue). Fur- 
thermore the plot suggests that Aj is the crucial parameter to 
trigger whether increasing interaction strength uq delocalizes 
or localizes the system. Technically, we evaluated systems of 
size L — 100 and averaged over 30 realizations. The dashed 
lines (black) sketch the two branches we present in Fig. 8 when 
additionally considering correlation. 

inverted and all even ones are left unchanged. The cor- 
responding operator in the basis reads 

= diag(...,-,+,-,-K, -,+,...) =W+ (27) 

and satisfies 1 = U^, i.e. U is a unitary transformation. 
In terms of annihilation (creation) operators one obtains 
Ucf ^ll+ = i-ycf'^ = e*'^Jc^+^ Now, the key observa- 
tion is that 

UHrnp{Uo)U+ = -H„,p{-Uo) (28) 

approximately holds when averaging over disorder real- 
izations (cf. Eq. (26)), since on average for each set {ej} 
there will be another one {—tj} when a = 0. More pre- 
cisely, if we fix the disorder potential {ej} the solution 
Hmp{-Uo) \E') = E' \E') is related to the solution of 
Hmp{Uo) \E) — E\E) with {— fj} by the identifications 
E' — —E and \E') =U\E) which obviously leaves phys- 
ical observables/localization measures invariant. 

Concerning the qualitative different localization prop- 
erties of the interacting bosons with respect to the range 
A/ we would like to provide some phenomenological ar- 
gument that may also draw a unified picture of delocal- 
ization by correlation and interaction. Starting from the 
case of short-range^'^ (Aj < 1), attractive interaction the 
enhancement of localization may result from the effect of 
lumping: Assuming sufficiently strong interaction com- 
pared to the kinetic energy {Uo > J), once the particles are 
nearby, they will be tightly bound in space and the quan- 
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FIG. 8: (color online) Impact of correlation (a > 0, o = 
as a reference) on the localization property of two interacting 
bosons. We picked out two representative branches (Aj = 0.04 
and Xj = 0.3) from Fig. 7 to demonstrate that, as expected, 
correlations tends to delocalize the system. Furthermore we 
note that the qualtitative result is again independent of the 
sign of the interaction and it seems that a feature from a = 
stays untouched: Namely, small interaction ranges A/ (weakly) 
delocalize with increasing interaction strength ug and with 
larger Aj states localize with respect to increasing uq. The 
numerical data were obtained at k = 5 from 30 realizations 
with L = 100 sites which we averaged over. 



turn dynamics should be governed by this localized be- 
havior. Due to the approximative invariance of the prob- 
lem to the interaction's sign the situation (surprisingly) 
stays unchanged even for repulsive interacting bosons. 
But according to Fig. 7 the effect of long-range interac- 
tion (Xj = 0.3) is qualitatively opposite to the short-range 
case and we numerically checked that the trend to de- 
localize with increasing interaction stays up to A/ = 1 
(not shown here). One may understand this behavior on 
the basis of the Fock space Hamilton matrix represen- 
tation Hmpjf (cf. App. A 2, Eq. (All) and Tab. I). Since 
disorder et and interaction Uij both contribute to the 
diagonal elements H^p jf we may effectively relate the 
resulting multi-particle Hamiltonian matrix to a corre- 
sponding one-dimensional, non-interacting system of size 
L{L + l)/2 with modified disorder ei'=i...z,(L-i-i)/2 and ex- 
tended dynamics beyond the nearest neighbor hopping^*. 
While short-range interaction just contributes to a few di- 
agonal entries of Hmpjf, an increasing interaction range 
Xj affects more and more of those matrix elements. Due 
to the deterministic character of the Uij (smooth, algebraic 
long-range decay) it will effectively Correlate the e^' of the 
non-interacting analoguous which we know to yield delo- 
calization for increasing long-range correlation. 

Finally, we want to include long-range correlation 
(a > 0) into our toy model of two interacting bosons 
to directly render the impact of correlation within our 
localization framework. In Fig. 8 we plot the uncorre- 
lated case from Fig. 7 for two different interaction ranges 
A/ as a reference and, therefore, it becomes clear that 



correlation among the onsite disorder yields delocaliza- 
tion — as we would expect from our previous experience 

so far (cf. Sect. II B as well as our phenomenological line of rea- 
soning from the preceeding paragraph). However, the feature 

that short-range [Xj = 0.04) interaction localizes and long- 
range (A/ = 0.3) interaction seems to delocalize the two 
bosons for increasing (rescaled) interaction strength uq is 
similar to the uncorrelated case. 

Again, the result seems to be independent of the sign of 
interaction. Since we have introduced correlation among 
the Ci our argument according to Eq. (28) should not be 
valid in general. Therefore we suppose the specific disor- 
der model to intrinsically fulfill the necessary assumption 
from above. Referring to the limiting case a — ^ oo we rec- 
ognize that ei becomes cos-like where — > — just shifts 
the onsite potential by a phase value tt. 

We would like to close our discussion by adding com- 
ments on the finite size analysis done as well as com- 
paring our investigations to similar literature, especially 
Ref. [42], who consider spinless fermions that interact 
when nearest neighbors to each other — an analogous 
to the (short-range) onsitc interaction scenario for bosons. 
Dukesz et al. employ the same disorder potential as us 
and study, among others, the interplay between long- 
range correlated disorder and interaction for two parti- 
cles and refer to it as the dilute limit. Their localization 
detection measure (NPC) is essentially the inverse of our 
IPR, but with the decisive distinction of being applied 
directly in Fock space — without any projection back to 
the lattice we advertised in Sect. HI A. Therefore we will 
be restricted to a comparison of qualitative features. 

Concerning the left upper panel of Fig. 6 in Dukesz 

et al. we notice (apart from minor deviation for small interac- 
tion strength denoted by A) the trend of enhanced localiza- 
tion (decreasing (NPC>) for increasing interaction strength. 
A feature that is confirmed by Fig. 7 for sufficiently 
small interaction range A/ (increasing IPR). Aside the non- 
interacting case {uo <-> A = 0) this qualtitative trend re- 
mains when long-range correlation is involved: Accord- 
ing to fixed a the different (NPC)-curves in the lower 
left panel of Fig. 6 decrease in magnitude. Moreover, 
each single curve supports enhanced delocalization for 
increasing correlation parameter < a < 4. 

To study effects of finite size it is actually wise to in- 
vestigate IPR , since our discussion on Eq. (6) suggests 
that it provides some notion on the effective number of 
occupied sites. It is perhaps that stage where one benefits 
from the proposed projection prescription of Sect. HI A 
again. While (NPC) measures occupation of the eigen- 

states \E) in Fock space our IPR ^ directly accounts for 
the projected probability distribution pi (cf. Eq. (21)) on 
the one-dimensional lattice with L sites. However, in 
order to establish a finite size analysis we depicted rep- 
resentative interaction and correlation parameter values 
(a, K, Uq, A/) and studied the IPR with increasing system 
size L in Fig. 9. For uncorrelated (a = 0.0) and corre- 
lated (a = 2.4) disorder we chose two different interaction 
strengths where there is no significant difference in the 
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FIG. 9: (color online) Finite size analysis of the IPR for representative model parameter values that specify the correlated 
disorder by a and the interaction by its strength uo and range A/. Due to the observed sign symmetry we restrict to the 
case of attractive interaction; the color encoding and point shape is kept as in Fig. 7. Instead of plotting the averaged inverse 
participation ratio IPR itself, it is more appropriate to consider its inverse since it directly relates to the number of occupied 
sites of the projected probability distribution pi (cf. Eq. (6) and Eq. (21)). 



IPR with respect to Fig. 8: a) uq = 1 and b) where the 
interaction range A/ significantly splits the IPR [uq = 9), 
respectively. 

In the case of uncorrelated disorder (left and middle left 
panel of Fig. 9) the averaged number of occupied sites 
L ~ IPR increases less than linear with linear increas- 
ing system size. Thus the extrapolation to the thermo- 
dynamic limit suggests localized states as for the non- 
interacting system^^. However, we would like to point 
out the increasing deviation of L between short-range and 
long-range interaction for strong interaction strength uq 
which we would like to phrase as enhanced delocalization 
by long-range interaction and which we suggest to keep 
in mind for further investigations of interacting particles 
in disordered media. If we turn to the strongly corre- 
lated regime {a > 2) we definitely observe a qualitative 

difference in the finite size analysis of L. Now the IPR 
seems to grow linearly in L and thus the relative num- 
ber of occupied sites L/L stays constant for L — > cx) 
assuming that the observed trend remains for L > 120. 
Hence strong correlation among the disorder potential 
delocalizes in analogy to the numerical experience from 
the non- interacting analysis, Sect. II B. But again, in the 
case of Uq = 9 the curves L{L) increasingly deviate for 
the scenario of short- and long-range interaction, which 
verifies our proposal of enhanced delocalization by long- 
range interaction and we scetched a potential (phenomcno- 
logical) reason by means of a (rough) correspondence to a 
non-interacting disordered system above. 

In summary, we hope that our data demonstrated that 
there is a complex interplay between correlation and in- 
teraction which by means of the Hamilton matrix struc- 
ture may arise from a similar origin, namely the correla- 
tion of the diagonal elements. Furthermore, both effects 
may compensate each other as seen for IPR|^_j^ Aj-oo4 
and IPR , „ , „ „ at un 7.5 in Fig. 8. The numeri- 

la— 1.6, A/— 0.3 ^ " 

cal study of the interacting bosons supports the indepen- 
dence of localization from the sign of interaction within 
our picture of multi-particle localization, and we hope 
that the rather academic treatment of two interacting 



bosons reveals in experiments with cold gases^'^; at least 
in principle in the limit of sufficiently low densities. Nev- 
erthless we want to remind that all ideas are based on 
rather phenomenological reasoning and numerical simu- 
lations where we are able to point out qualitative trends 
only. Therefore our results can just pave the way for a 
more profound understanding on theoretical grounds. 



IV. CONCLUSION & PERSPECTIVES 

To conclude our investigations we briefly summarize what 
was achieved within our study of the one-dimensional 
Anderson model. In the first place we analyzed two 
mechanisms that yield delocalization, namely a corre- 
lated disorder potential and mutual interaction between 
two bosons. By reviewing a known model of correlated 
disorder from the perspective of exact numerical diago- 
nalization and different (well-established) localization mea- 
sures, we showed that the process of delocalization is 
much more complex than observed up to now: Namely, 
there are reasonable arguments to consider a = 1 as im- 
portant for the process of delocalization by correlation. 
Nevertheless the general trend that increasing correlation 
yields delocalization was confirmed and localized states 
approximately vanish for a > 2. In addition we were able 
to numerically establish the full parameter dependence of 
the model system by means of the nearest neighbor distri- 
bution (NND) and therefore we went beyond the analysis 
that focuses on states in the band center. The results sug- 
gest that the system eventually becomes delocalized for 
sufficient large correlation, independent from the disorder 
strength k which is in accordance with Bloch's theorem. 
Nevertheless, we extensively discussed the usage of the 
NND and tried to convince that it is far less obvious to 
utilize it as an localization detection measure. 

Therefore we introduced a general idea to define lo- 
calization for multi-particle states of bosons by means of 
the phenomenology of experiments with cold quantum 
gases. The examination of two long-range interacting 
bosons on a finite, one-dimensional lattice with (perfect) 
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FIG. 10: Illustration of autocorrelation (cf. Eq. (A3)) among the onsite disorder e^, Eq. (3). We evaluated 5000 points (ei, ei+i) 
with A'^ — 10^ (cf. Eq. (3)). An increasing correlation parameter a deforms the rotationally invariant density of points (ei, ej+i) 
to a linear dependence between and ei+i which can be quantitatively detected by the (linear) autocorrelation C{d). Rotational 
invariance yields C(l) = and linear dependence between and e^+i implies C(l) = ±1. Note, that even in the case a = 
our model disorder is not totally homogeneous in [—0.5, 0.5] x [—0.5, 0.5]. 



disorder confirmed the conjecture that the localization 
property does not dependent on the sign of interaction. 
We showed that this phenomenon leads back to a dis- 
crete symmetry arising when averaging over disorder re- 
alizations. Finally the complex interplay of correlation 
and interaction was studied numerically and by means of 
phenomenological reasoning on the basis of the Hamilton 
matrix structure we argued how to relate the impact of 
interaction to delocalization by correlation known from 
the non-interacting system. 

Although our contribution provides some new insights 
to the phenomenon of localization, unsolved aspects re- 
main. More precisely, a solid theoretical description of 
the crossover from localization to delocalization when 
tuning the (nearly) perfect disorder to the Bloch-like situ- 
ation of a highly correlated potential is desired. Moreover 
a more profound understanding of the numerically ob- 
served effects of interaction on localization is imperative 
to obtain further insight to delocalization/localization in 
the presence of many interacting particles. 
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Appendix A: Numerical Details 

This appendix is dedicated to computational details of 
our work. We especially address the correlated disorder 
used in the main part and an explicit calculation of the 



Hamiltonian matrix of two interacting bosons for an ar- 
bitrary (finite) number of lattice sites. 



1. Correlated disorder 

Here, we provide some properties of the correlated disor- 
der potential and as a general remark we mention that 
all our numerical results were averaged over a number of 
disorder realizations. If we exactly follow the prescription 
Eq. (3) to obtain L values with correlation parameter 
a the quantity 

6 = max ei — min , i ~ 1 . . . L (Al) 

i i 

will not be constant in general. To ensure the same dis- 
order strength k for all disorder realizations and system 
sizes, we always normalize a given chain of values by 

A A 

Ci ^ Ci- mme^ -> —e, -> - — (A2) 

i d I 

to have ti £ [— y, ^]. See also Refs. [43,44] for a related 
discussion. 

Since we mentioned that a = refers to nearly per- 
fect disorder, we provide an intuition on that statement 
in Fig. 10 where we illustrate increasing correlation by 
plotting ei versus e^+i for different (increasing) a-values: 
For a — the point density is approximately rotationally 
invariant around (ei,ei_|_i) = (0,0) and for the autocor- 
relation function 

C{d) = hm ^ (A3) 

where 

{^^)L ^ i E ^ (e?)^ - , (A4) 

i=l 
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we find^i C(l) « for a = 0. By means of Eq. (A3) 
we can state that the disorder is uncorrelated. But 
as obvious from Fig. 10 the whole plane (ei,ei+i) £ 
[-0.5,0.5] X [-0.5,0.5] (A = 1) is not uniformly filled 
with points (ei,ei-(-i) and therefore an arbitrary onsite 
potential value can not obviously be followed by any 
possible e^^i. Thus we refer to the case a = as weakly 
correlated. On increasing correlation we observe that the 
(ei,ei+i) become more and more linearly correlated, i.e. 
Ci ~ const. ■ Ei+i and thus C(l) « 1. 

Let us finally sketch one feature of the that arises 
when approximately computing C(d) to demonstrate the 
long-range character of the correlated disorder. Let 

/•OO 

e{xi) / dk k^"^^ cos[kxi + Xi ^ iAx (A5) 

Jo 

be the continuous version of Eq. (3) and 

C{d) cx lim 4 / dx e{x)e{x + d) (A6) 

X-foo X Jq 

that of Eq. (A3). Hence we obtain 

1 r°° 

C{d) (X - dk fc-" cos(fcd) (A7) 

2 Jo 

which can be analytically solved for a E (0, 1) and we 
end up with 

C{d) oc < a < 1 . (A8) 

Hence we have a direct argument at hand why the ej are 
long-range correlated: Its spatial autocorrelation decays 

algebraically (at least in the continuum limit for a e (0, 1)) SUch 

that increasing a yields a slower decay with distance d. 
Within this picture a = 1 is also highly correlated and 
one may argue that delocalization occurs below a — 2 
which is qualitatively supported by our localization mea- 
sures although the LL agrees pretty well with other re- 
sults from the literature (cf. discussion in Sect. 11 B). 



2. Matrix representation of two interacting bosons 
in the modified/interacting Anderson model 

For the numerical results on two interacting particles we 
developed a rather general scheme of the matrix represen- 
tation in Fock space which allows us to set the numerics 
for an arbitrary number of sites. Moreover the matrix 
is banded with a width ~ L which helps to increase the 
efficiency of the diagonalization routine. As we were in- 
spired by an experimental setup, we used hard-wall/open 
boundary conditions, i.e. there is no (kinetic) hopping el- 
ement J (cf. beginning of Sect. (I)) from sitc i = 1 to site L 

and vice versa. 

We obtained the matrix structure presented in Tab. I 
where we labeled and ordered the Fock states according 
to 





[200., 


..0) 


|2)^ = 


[110., 


..0) 


|3)^ = 


[101., 


..0) 


\L)p = 


[100., 


..1) 


1)^ = 


[020.. 


..0) 



|L(i + l)/2)j, = [000... 2) (AlO) 

with \n1n2 . ■ ■ ni) indicating the number Ui of bosons at 
site i and the matrix structure was exemplarily evalu- 
ated for L — Q. Entries x denote (constant) kinetic off- 
diagonal values and x refers to the diagonal elements that 
are determined by the correlated disorder terms Ci and 
the interaction potential Uij. If we write the Hamiltonian 
matrix elemtens as 

Hmp,ff' = {f\H„ip\f')p (All) 

the generalization to an arbitrary number of sites L reads: 



□ diagonal element H„ip J f contributions 

f: 1 2 ... L L+1 L+2 ... L(L+l)/2-l L(L+l)/2 

QC+c,;: ei + ei Ci + £2 ■•• ei + el £2 + £2 £2 + ^3 ■•■ <^L-i + <^L <^L + 

U,,,ctc+c,cf UUO) C/f(l) ... C/±(i-l) C/±(0) [/±(1) ... C/±(l) [/±(0) 

□ first off-diagonal elements Hmpj f+i/ J 

Jct^j^cf. V2 1 ... V2 1 ... 

□ remaining off-diagonal elements H„ipj f 

Jc^j^q: all remaining off-diagonal elements have magnitude J. In the jth off-diagonal (j > 1) there are j 
non-zero elements: Hmpjf+j with 

f^L~j:l^,r^L-j{j + l)/2, ... , L-j(j + l)/2-(j-l) . 
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(A9) 



TABLE I: Matrix structure of the multi-particle Hamiltoriian Hmp in the 2-particle basis, Eq. (AlO) with L = 6 sites. 
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As usual in condensed matter physics, the flow is 
parametrized by the system size. 

In our numerics with L ~ 10'^ — 10^ we choose No = 10^ 
and ensure that our results are not altered for > A'^o by 
means of finite size analysis, i.e. we do not monitor any 
qualitative deviation of the data for A'^ > Nq. 
''^ On the basis of coupling two harmonic oscillators the no- 
tion of level repulsion becomes a bit more explicit. In fact 
Hs from Eq. (1) just mimics such a system with ground 
state energy difference Ae = |ei — £2] and (ground state) 
coupling energy J between the two oscillators labeled by 
i = 1,2 in the most simplest case of L — 2. By inves- 
tigating the quantum physics of a single particle in such 
a setup with ei = one observes that the energy differ- 
ence of the system's eigenstates reads AE = — _E2| = 
2\/Ae2/4 + J2, i.e. in the uncoupled case J = the eigen- 
states \Ei =0) and \E2 = Ae) are localized in one of the 
harmonic oscillators and by tuning the external parame- 
ter Ae — >■ their corresponding eigenenergies Ei become 
degenerate. By means of A_E this is impossible for J 7^ 
where both eigenstates , I-E2) are distributed among 
the two coupled oscillators. 

Despite similar in notion the Lyapunov exponent 7 in local- 
ization theory has not to be confused with its counterpart 
in nonlinear dynamics*^ 7 (where it is often denoted as A). 
Here, 7 is equivalent to the exponential decay rate of -tp^ 
(cf. Eq. (4)) for sufficient large systems. 
In particular J^i^ dE DOS(£') is the number of states \E) 
with El < E < E2 and rescaling the DOS simultaneously 



with the spectrum {E} is understood to fulfill the normal- 
ization condition 1 = / dEDOS{E). 

As known from standard text books, e.g. Ref. [4] or F. 
Dyson's work Ref. [30] one assumes a) the joint probabil- 
ity p(//rand) of the raudom matrix ensemble -ffrand being 
invariant under orthogonal transformations and b) one re- 
quires the statistical independence of all Hamilton matrix 
elements which does not hold in our case of i^rand Hs '■ 
Despite the random potential values the kinetic hopping 
energy J is deterministic/constant. 

According to their notation k translates to W and a to p. 
While we directly rely on Eq. (3) to produce correlated 
disorder (see also App. Al), Ref. [11] performs some sort 
of Fourier filtering method to obtain long-range correlated 
Ei-values. 

We choose the interaction Uij to dependent on the system 
size L to ensure the Constance of the interaction within 
the picture of the continuum limit described above when 
L tends to infinity. 
^® Constant filling means n L and together with 

dimZ/mp = (n + L— iy./{n\{L— from elementary com- 
binatorics and Stirling's approximation logs! ~ a; log a; — a; 
we obtain dim_H"mp ~ for L ^ 1 and a ~ 1. 
We would like to note that there is a minor inaccuracy 
in notion here. Although the algebraic decay of the inter- 
action potetial Uij for |i — j| ~ L is always long-ranged 
compared to e.g. an exponential decrease, in the following 
the terms short- and long-range will refer to A/ ^ 1 and 
A/ ~ 1, respectively. 

We would like to remind that Anderson's initiating paper^ 
did not restrict to the case of nearest neighbor hopping. 
That was also the reason why we mentioned to deal with 
a certain subclass of Hamiltonians, Eq. (1), in Sect. I. 
^® This conclusion is opposite to Dukesz et al. who present 
their results in the lower right panel of Fig. 7 in Ref. [42] . It 
would be interesting to extend there study to larger system 
sizes beyond L = 60. On the other hand the deviation of 
the results suggest that a proper projection to the lattice 
is perhaps crucial in order to investigate localization of 
multiple particles. 

For this purpose it is obviously necessary to produce cor- 
related long-range disorder with algebraic power spectrum 
which could be a challenging task for experiments. Via Fes- 
hbach resonances'^^' it is possible to tune the interaction 
strength uq, but the variation of the range A/ is perhaps 
difficult to establish. 

In general, correlation between et and €i+d is given by C{d). 



